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ABSTRACT 

We investigate the high-energy (HE) (< IGeV) emission from the pulsar PSR B1509- 
58 and its relation to the radio emission in the 1.4GHz frequency band. The role of 
the quasi-linear diffusion in producing the pulsed HE radiation is investigated. We 
show that by means of the cyclotron instability the relatively low frequency waves 
excite, which due to the diffusion process influence the particle distribution function 
and switch on the synchrotron emission mechanism. We argue that the coincidence of 
HE main peak and the radio pulse is a direct consequence of the fact that the high 
and low frequency radiation is produced simultaneously in a local area of the pulsar 
magnetosphere. In the paper we also consider the absence of the radio counter pulse 
and explain this fact. 
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1 INTRODUCTION 

Recent development in tlie field of high energy 7-ray as- 
tronomy stimulates new theoretical approaches and meth- 
ods in modern astrophysics. Last decade the bulk of ob- 
servational results in the high and very high energy do- 
mains comes from several telescopes: High Energy Stereo- 
scopic System (HESS), Very Energetic Radiation Imaging 
Telescope Array System (VERITAS), Major Atmospheric 
Gamma-ray Imaging Chcrenkov (MAGIC) Telescope etc. 
These instruments have detected several sources at GeV 
and TeV energies le.g. Mariotti"1l2010l : lAbramowski et al. I 
I2OIII : lOng fc Mariotti II2OIOI ') and Fermi Large Area Tele- 
scope (Fermi-LAT) in turn, an nounced the first c atalogues 
of ac tive galactic nuclei (AGN) (|Abdo et al. l20l"0al ) and pul- 
sars (|Abdo et al. I|201l l). respectively. 

In the context of the present paper a special interest 
deserves the detection of the HE < IGeV 7-ray pulsed emis- 
sion from the young and energetic pulsar PSR B1509-58 by 
the Fermi-LAT with a light curve presenting two symmet- 
ric peaks separated by phase 0.37 ± 0.02 (see Fig.l from 
lAbdo et al. I (|2010bll '). The 7-ray signals are offset from the 
main radio peak observed at L4GHz by phases 0.96 ± 0.01 
and 0.33 ± 0.02. The first 7-ray peak is almost coincident 
with the radio pulse. Particularly, leading the main radio 
pulse by phase 0.04 ± 0.01. The results of Fermi-LAT are 
obtained for the energy interval 30MeV-100GeV and it is 
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shown that up to IGeV the emission reveals pulsations. We 
focus on explaining the pulsed HE 7-ray emission of this 
source in IMeV-lGeV interval. 

According to the observational evidence, the main 
pulses in the radio (1.4GHz) and HE (< IGeV) bands 
are almost coincident, meaning that they are generated 
simultaneously, in a certain localized area of the pulsar 
magnetosphere. In this context the quasi-linear diffusion 
(QLD) could be an eflicient mechanism providing the 
coincidence of pulse s in the aforementioned domains. 
iKazbegi et al. I (| 19921 ) showed that if the magnetic en- 
ergy density exceeds that of plasmas the low frequency 
cyclotron modes are excited. In turn, the unstable cy- 
clotron waves interact with resonant particles by means 
of the QLD and as a result they acquire the pitch 
angles, leading to the effici ent synchrotron emission pro- 
cess as i t was argued b y iLominad ze et alTI (| 19791 ) and 
Machab eU fc Usovl (|l979l i For high Lorentz-factors the 
synchrotron emission comes in the HE domain. In a series 
of papers, the mechanism of QLD was applied to active 
galact ic nuclei (|Osmanov fc Machabeli I I2OI0I: lOsmanov I 
2010allbl. I2OIII') and pu ls ars jMalov fc Machabeli ' '2OOII: 
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20091 . 12OIOI : IChkheidze et al. IbOllT) . respectively Recently 



we have applied this mechanism to explain the pulsed 
very high energy (> 25GeV) ra di ation o f the Crab pulsar 
(IMachabe li fc Osmanov ] I2OO9', '2OIOI; IChkheidze et al. I 
[2011) and it's pulse-ph ase coincidence with the optical 
signals jAliu et al. Il2008l h 
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In the present paper we implement the mechanism of 
the QLD to the pulsar PSR B1509-58. The observations 
confirm that the signals of the HE 7-rays (< IGeV) are 
almost coincident in phase with the radio signals (1.4GHz). 
Thus, we argue that they are produced simultaneously in 
one location of the pulsar magnetosphere. 

The paper is organized as follows. In Section 2 we con- 
sider our model of the quasi-linear diffusion and apply it to 
the pulsar PSR B1509-58, explaining major characteristics 
of emission; in Sect. 3 we discuss the model and in Sect. 4 
we summarize our results. 



2 MODEL 

According to the standard model of pulsars, it is be- 
lieved that the pulsar's magnetospheres is filled by electron- 
positron plasma, which consists of three components. The 
so-called primary beam particles, which uproot from the 
neutron star's surface and are described by the Goldreich- 
Julian density l|Goldreich fc Julian 111969^ in due course of 
time accelerate in the gap electric field and reach a cor- 
responding thres hold when the cascade pair crea tion be- 
comes inevitable (jSturrock Ill97ll : iTademaru Ill973h . There- 
fore, the pulsar's magnetosphere is filled by the primary 
beam (7 ~ 75), the bulk of plasma (7 ~ 7,,) and the s o- 
called tail particles (7 ~ 7t) (see Figl. from I Arons I l|l98ll) ). 
The distribution function of the electron-positron plasma is 
one-dimensional, as at the pulsar surface any transverse mo- 
menta of relativistic electrons are lost in a very short time 
10~^''s) via synchrotron emission in very strong magnetic 

fields. 

According to the work of iKazbegi et al. I l|l992l ). by 
means of the anomalous Doppler effect plasmas in the pulsar 
magnetosphere undergo the cyclotron instability. In partic- 
ular, due to the curvature radius of the field lines, p, the 
beam particles are forced to drift perpendicularly to plane 
of the curved field line with the following velocity 



pUJB 



(1) 



where c is the speed of light, is particle velocity along 
the field line, lob = eB/mc is the cyclotron frequency, B is 
the magnetic induction, e and m are the electron's charge 
and the rest mass, respectively. The distribution function 
is one-dimensional and anisotropic and plasma becomes un- 
stable. Both of these factors (the one-dimensionality of the 
distribution function and the drift of particles) might cause 
generation of eigen modes in the electron-positron plasma if 
the following resonance condition is satisfied 



uj — k„ V„ — kxUx — — — 

II II 



(2) 



Here, fc|| is the longitudinal component of the wave vec- 
tor and kx is the wave vector's component along the drift. 
From Eq. ((2]) one can see that the pure transversal modes 

i t-waves) characterize d by the following dispersion relation 
Kazbegi etHII 19921 ). 



u!t ~ kc {1 — S) . 
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(3) 



A-KUpe'^ /m is the plasma frequency and Up is the plasma 
density. In the framework of the present model, the afore- 
mentioned waves, by means of the QLD, influence the reso- 
nant particles leading to the diffusion of particles along and 
across the magnetic field lines. Consequently, particles ac- 
quire pitch angles and start to radiate in the synchrotron 
regime. 

It should be mentioned that we study the process of 
the QLD close to the light cylinder surface, where, as it 
was investigated, effe cts of centrifugal accelerat ion are very 
efficient for pulsars (jOsmanov fc Riegerl [2009') and AGN 
jOsmanov et al. II2007I ). Particularly, for explaining high and 
very high energy emission of pulsars one needs a mecha- 
nism providing very high energies of electrons, which is not 
pr ovided by the gap accelerat ion process. As it was shown 
in lOsmanov fc Rieger I l|2009l ) , the centrifugal acceleration 
process is quite efficient in the light cylinder zone and can 
provide the reacceleration of the primary beam particles to 
reach the Lorentz factors of the order of 10^. This in turn 
leads to the HE emission patt ern, that is in a good agree ment 
with the observational data (jOsmanov fc Rieger IlioO^ . 

In order to describe the process of the QLD, one 
has to take into account as the dissipative factors, which 
try to decrease the pitch angles, as the diffusion attempt- 
ing to increase them. When relativistic electrons emit in 
the synchrotron regime, they undergo the so-called reac- 
tion force, having as longi tudinal as transversal components 
l|Landau fc Lifshitz |[l97ll ) 



Fi_ = 



,Fii = 



2 



(4) 



where a = 2e^aj^/3c^, px and p\\ are transversal and longi- 
tudinal components of momentum, respectively. 

On the other hand, the diffusion process, under cer- 
tain conditions, might balance the corresponding dissipa- 
tion, leading to the quasi-stationary state (djdt — 0), which, 
as is shown by Chkheidze et at. (2011), for d/dp± » d/dp\\ 
results in the following kinetic equation governing the afore- 
mentioned mechanism 



7^(pxF,/) = /- 
dpi_ dp± 



P±Dj 



df_ 

dp± 



where / is the distribution function of electrons. 



8c 



(5) 



(6) 



is the diffusion coefficient and Ek is the electric field. In or- 
der to estimate the value of \Ek\^ , we assume that half of 
the beam energy density, mc'^nt"ft/2 [nt is the beam den- 
sity), converts to the energy density of the waves |_Efcpfc. 
Therefore, an expression of l-E^p writes as follows 



\Ek. 



mc Ubjb 
2^ ■ 



The solution of Eq. @ is given by 



f{p±) = Cexp 



Fi_ 



dp± = Ce 



where k denotes the modulus of the wave vector, cup = 



where 

P±o = 



(7) 



(8) 



(9) 
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and for the mean value of the pitch angle we obtain 



(10) 



2.1 Radio emission 

According to the present model, the generation of the syn- 
chrotron emission is provided by the feedback of the low 
frequency cyclotron waves on the resonant particles. It is 
easy to show from Eqs. ((2l [3| that the unstable cyclotron 
wave s will have the following frequency (|Machabeli fc Usov I 
[l92i) 



X IQ-^'GHz. 



(11) 



The estimations show that for reasonable magnetospheric 
parameters (we assumed that 7p ~ 1 and 7;, ~ lO'^ 
l|Osmanov fc Rieger I [2OO9I )) the anomalous Doppler effect 
leads to the excitation of the radio waves (1.4GHz), which in 
turn influences the particles and switches on the synchrotron 
ra diation proce ss giving the HE 7-ray emission. According 
to lAbdo et al. I (|2010bi ) it is shown that the main pulses of 
the HE (< IGeV) 7-rays and the radio emission come with 
the same phases. We argue that this observation feature is 
caused by the fact that the emission in both domains are 
produced simultaneously and in one location of the magne- 
tosphere - close to the light cylinder surface. 



2.2 Spectrum of the high energy emission 

In this subsection we are going to explain the observed 
spectrum of the HE 7-rays. Following the method devel- 
oped by Ginzburg (1981), one can show that the flux of the 
synchrotron emission is given by the following expression 
(Chkheidze et al. II2OIII ) 



oc 



Kr^l-i{z)dz 



dp„ , (12) 



where /^^ (p|| ) is the longitudinal distribution function of elec- 
trons. The major difference from the model of Ginzburg is 
that according to our scenario the synchrotron emission is 
generated in the outer part of the pulsar magnetosphere via 
the QLD, which inevitably restricts the values of the pitch 
angles. On the other hand, according to Ginzburg (1981), 
along the line of sight the magnetic field is chaotic leading 
to the broad interval of the latter i/> € (0;7r). This differ- 
ence is principle between the standard approach and the 
method presented in this paper. Therefore, in Eq. p2p in- 
stead of tp we put it's average value (see Eq. (|10|l 'l. which, 
unlike the standard model, is the dynamical parameter and 
participates in integration. 

In order to find out the spectrum of the HE emission, 
one needs to know the behaviour of the distribution function, 
/ll(Pll). For solving this problem, we consider the equation 
governing the evolution of / (p,,) l|Chkheidze et al71l201ll ) 



dt 



d 

dp, 



i2„2 7^1/2 



2 J- 
P±of„ 



(13) 



We are interested in the stationary state {d/dt — 0), of the 
QLD, which reduces the aforementioned equation and leads 
to the following form of the distribution function 



II 



(14) 



On the other hand, the cyclotron noise is described by 
the equation 



d\Ek\ 
dt 



(15) 



where the incremen t of the instability, F, is given by 
l|Kazbegi et al. II 19921 ) 



(16) 



Here lui, = ^ Aimbe'^/m is the plasma frequency of the res- 
onant (beam) electrons. 

By taking into account Eqs. (I15 | 16[) one can re duce Eq. 
(fT3)) to the following form l|Chkheidze et al. Il201lh 



(17) 




where 



3 TT^C- 



(18) 



and /||^ is the initial distribution function of particles. Since 

the electron number density behaves with distance as 1/r'', 

it is clear that one has the following condition, f ^ f , 

° II l|o' 

that reduces Eq. ^\ 




+ f =0. 



(19) 



As we can see the function Ek{p\\) drastically depends 
on the form of the initial distribution of the primary beam 
electr ons. According to the work, (jGoldreich &: Julian I 
I1969D . a spinning magnetized neutron star generates an elec- 
tric field which extracts electrons from the star's surface and 
accelerates them to form a low-density (nj, — B/Pce) and 
energetic primary beam, which in turn, is reaccelerated in 
the light cylinde r zone and might reach qu ite high Lorentz 
factors, ~ 10^, (|o smanov fc Rieger I [20091) . We only know 
the scenario of creation of the primary beam, but nothing 
can be told about its distribution, which drastically depends 
on the neutron star's surface properties and temperature. To 
our knowledge there is no convincing theory which would 
predict the initial form of the distribution function of the 
beam electrons. Thus, we can only assume that the beam 
electrons have the power-law distribution, f oc p~", and 

for l-Efcp we will obtain the following behaviour 



I TTi 1 2 3 — 2n 

\Ek\ ccp,, 



(20) 



In order to simplify Eq. (I12|) . we replace the integration vari- 
able by X = e/eccv — 3-6 X 10^*i^7~^'^'', and taking into 
account Eqs. (|6I9I10I14I20P the expression of the fiux will be 



R (X e~ 



xi-- 



Kr,/z{z)dz 



dx. 



(21) 
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Figure 1. Spectral energy distribution of PSR B1509-58 (dashed- 
line corresponds to the model spectrum e~^'*^ea;p[— (e/0.059)^''^] 
and the solid- line to e~^'*'^e3;p[— (e/0.081)] obtained by 
iPilia et a. 1 1 I2OICII ) from a fit of pulsed fluxes from soft to hard 7- 
rays. The round points represent COMPTEL and square p oints 
AGILE observations (Kuiper et al. 1999; Pil ia et a. II2OICII). re- 
spectiv ely. The arrow represents Fermi upper limit lAbdo et al. I 
l2010bh . 



As we see from this equation the HE spectral index is given 
by the foll o wing expression (2 — n)/(4 — n). According to 
iPilia et a. I (|2010l ) the observed HE pulsed emission of PSR 
B1509-58 in the energy domain (0.001 - l)GeV is best de- 
scribed by a power-law plus cutoff with the spectral index 
equal to 1.87. When n = 6.3 the spectral index of the syn- 
chrotron emission equals 1.87 (see Eq. (|21[) V If we as- 
sume that the energy of the emitting electrons vary between 
'jrnin — S-IO'' and 7maa; — 3.5-10'', oue cau show that the in- 
tegral (|21|) can be approximately expressed by the following 
function 



(0.059) 



(22) 



As we can see our emission scenario predicts the exponential 
cutoff, with the cutoff energy 59MeV (see Fig. [T]). 



3 DISCUSSION 

It is clear that the present model provides the self con- 
sistent explanation of several interesting observational fea- 
tures of the pulsar PSR B1509-58. In particular, we have 
already seen that for reasonable magnetospheric parame- 
ters the QLD can support the HE pulsed emission up to 
IGeV. Acco r ding to the observations, as it was shown by 
lAbdo et al. I (|2010al '). the main pulses of the HE domain are 
almost coincident with those of the radio emission on the fre- 
quency 1.4GHz. In the framework of our model this might 
mean that radiation in both domains are produced almost 
simultaneously and in one location of the pulsar magneto- 
sphere. We have shown that by means of the cyclotron in- 
stability, the radio waves are excited, which influence the 
resonant particles and cause their diffusion as along also 
across the magnetic field lines. Consequently, they acquire 
the pitch angles and the synchrotron emission process is 
switched on. This is the reason, why the HE 7-rays and the 
low frequency waves are strongly connected. On the other 
hand we neglect the role of the curvature radiation, because 




Figure 2. Emission geometry of PSR B1509-58. 



otherwise the main pulse of the HE emission would be wider 
than it is observed. This automatically means that the cur- 
vature radius of the magnetic field lines are large enough 
to insure the negligible role of the curvature radiation pro- 
cess. As it was shown by Osmanov et al. (2008,2009) such a 
configuration of the field lines might be provided by means 
of the so-called curvature drift instability, that makes the 
magnetic field lines rectify efficiently in the very vicinity 
of the hght cylinder surface. In particular, as it was shown 
in the aforementioned articles, the curvature drift (see Eq. 
IT])) will create the corresponding current, leading to the 
creation of additional magnetic field. On the other hand, 
this toroidal field, by means of the centrifugal effects, will 
efficiently amplify, changing the overall configuration of the 
magnetosphere. It has been found that the curvature drift 
instability has two general modes, one of which makes the 
field lines more curved and the other makes them rectify. 
We assume that since the HE profile of the main signal is 
relatively narrow, the curvature radiation should be negli- 
gible and therefore, the field lines must be almost straight, 
which is provided by the curvature drift instability. 

In the framework of the present paper we provide the- 
oretical confirmation of the measured HE spectrum of PSR 
B1509-58 in the energy domain (0.001 - l)GeV. In partic- 
ular, our model predicts the power-law plus an exponen- 
tial cutoff, which provides an explan ation of the observe d 
decrease of the flux above lOMeV (|Abdo et al. I l2010bl ). 
If the Lorentz factor of the resonant beam electrons are 
7i, — 3 • lO'' — 3.5 • 10^, then our theoretical spectrum 
Fe oc e~^'*^ea::p[— (e/0.059)'^'^]. As we can see the model au- 
tomatically yields the exponential cutoff, with the cutoff en- 
ergy 59MeV. 

Another important feature of the high energy 7-rays is 
that unlike the radio emission, apart from the main pulse 
it also has the second pulse, separated from it by the phase 
A(^ = 0.37 ±0.02. 

The observationally evident fact that the HE emission 
reveals two peaks might be a direct result of an emission 
geometry. In particular, we detect the pulsar emission twice 
per one rotation if the HE emission cone is almost perpen- 
dicularly inclined with respect to the rotation axis. Accord- 
ing to the observations the main 7-ray pulse leads the radio 
pulse by phase 0.04, which indicates that the emission cones 
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of the HE and radio domains are not spacially coincident 
(see Fig. [2]). In Fig. [2] cone 1 represents the HE emission 
cone, that has a common area (area 1 and area 2) with the 
radio emission cone (cone 2). Then it is clear that if the 
emission comes from area 1 which is parallel to the line of 
sight, we will see the emission profiles for the radio and HE 
bands, respectively. Due to star rotation and the correspond- 
ing emission channels, the counter HE radiation cone will 
have a position convenient for detecting the counter signal, 
whereas the counter radio cone will be significantly shifted, 
and the line of sight will be out of it. Therefore, as a result, 
we see the signal and the counter signal of HE emission and 
do not see the counter radio peak. 

This approach explains another feature of the HE emis- 
sion as well. As it is evident from the observations, the 
counter peak is wider with respect to the main pulse. In par- 
ticular, from Fig. (2] it is clear that the time during which 
the line of sight stays inside the cone 1 is shorter than the 
time scale it stays inside counter HE cone. 
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